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Control of negative thermal expansion is a fundamentally interesting topic in the negative thermal 
expansion materials in order for the future applications. However, it is a challenge to control the negative 
thermal expansion in individual pure materials over a large scale. Here, we report an effective way to control 
the coefficient of thermal expansion from a giant negative to a near zero thermal expansion by means of 
adjusting the spontaneous volume ferroelectrostriction (SVFS) in the system of PbTi0 3 -(Bi,La)Fe0 3 
ferroelectrics. The adjustable range of thermal expansion contains most negative thermal expansion 
materials. The abnormal property of negative or zero thermal expansion previously observed in 
ferroelectrics is well understood according to the present new concept of spontaneous volume 
ferroelectrostriction. The present studies could be useful to control of thermal expansion of ferroelectrics, 
and could be extended to multiferroic materials whose properties of both ferroelectricity and magnetism are 
coupled with thermal expansion. 

The thermal expansion compatibility of different component is a one of key problems for modern devices, 
such as thin film, multilayer chip capacitors (MLCCs), solid oxide fuel cell (SOFC), thermoelectric materials, 
and high temperature piezoelectrics 14 . Thermal stress could give rise to material failure due to undesirable 
mismatch of coefficient of thermal expansion (CTE) of different components. The control of thermal expansion is 
an important topic for the material design. Negative thermal expansion (NTE) was recently found in some 
materials in the past two decades 5 " 11 . NTE materials have a promising potential application on tailoring CTE 7 . 
The overall thermal expansion can be tailored by either formation of composites or chemical modification of 
single-phase materials. For example, some composites were fabricated, such as ZrW 2 0 8 with common materi- 
als 1213 . The progress is, however, hampered mainly due to the poor thermal stability of ZrW 2 O g . It decomposes at 
a relatively low temperature (777°C) 7 . If the control of CTE could be achieved in single-phase materials, draw- 
backs in composites such as low binding capacity, thermal stress, or chemical reaction at interface, could be 
effectively overcame. 

The control of NTE is an interesting topic for the studies on NTE especially of individual pure materials. Many 
kinds of chemical modifications have been carried out on various NTE materials. For example, ZrW 2 0 8 based 
single-phase materials show little dependence of NTE on compositions and behave a very small varied range of 

CTE (a v = -2.2 2.6 X KT^C" 1 ) by chemical modifications 14 . The control of NTE of ZrW 2 O s is restricted 

by the fact of few suitable substitutions and low solubility limit. A single-phase ceramic material of 
Al 2x (Hf,Mg) 1 . x (W0 4 ) 3 was found to have a bulk CTE range of -0.23 ~ 0.45 X Kr 50 C _1 , where a zero thermal 
expansion (ZTE) was also achieved 8 . Orthophosphate group of NaZr 2 (P0 4 ) 3 (NZP) of fast ionic conductivity 
have been widely studied with various substitutions at different sites to exhibit a CTE range of —0.17 ~ 0.19 X 
10~ 5O C~ 115 . The effect of cation size could play an important role in the thermal expansion property of NZP-based 
materials 16 . More recently, anti-perovskite Mn 3 AN (A = Zn and Ga) nitrides were found to have a giant NTE 
which was suggested to correlate with the coupling between the T 5g antiferromagnetic spin configuration and 
lattice 1718 . Many kinds of substitutions have been tried at three crystallographic sites in order for adjusting 
thermal expansion of Mn 3 AN. ZTE could be established in, for example, nanocrystalline Mn 3 Cuo. 5 Ge 0 .5N below 
— 55 °C 19 , and Mn 3 (Ga 0 .5Ge 0 .4Mn 0 .i)(Ni_ x CJ in a temperature range of —83 ~ 0 °C 20 . The further studies on the 
complex interplay between magnetic exchange, cation radius, and lattice would allow more precisely control NTE 
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of Mn 3 AN materials, especially for a wide temperature range. The 
challenge still remains for the control of single-phase NTE materials. 
Especially, the control of NTE over a large range would benefit the 
practical applications of NTE materials. 

Among the family of NTE materials, there is one branch of ferro- 
electrics which behave a general phenomenon that unit cell volume 
contracts during ferroelectric-to-paraelectric phase transition, 
such as perovskites of PbTi0 3 (PT) 21 " 23 , BaTi0 3 24 , BiFe0 3 25 , and 
(Bi,La)Ni0 3 9 , tungsten bronze PbNb 2 0 6 26 , and semiconductor 
GeTe 27 . It is important to note that all ferroelectrics except PT (RT 
~ 490°C) 23 exhibit NTE property within a narrow temperature 
range, e.g. 50 °C for BiFe0 3 25 . Such large NTE temperature range of 
PT offers us an opportunity to control the NTE property. Our pre- 
vious studies have revealed that the CTE can be adjusted by chemical 
substitutions for either Pb or Ti cations 21 " 23 . The temperature 
dependence of spontaneous polarization (P s ) and lattice dynamics 
have approved that the ferroelectric behaviour plays an important 
role in the NTE nature of PT-based materials 21,22 . For future applica- 
tion of the NTE PT-based materials, however, one of main problems 
still remains. That is how to easily control the NTE over a large scale. 
Here, we report an effective way to control the NTE of ferroelectrics 
over a giant scale through a simple method by adjusting the ratio of 
cations with different ferroelectric activity which is defined as the P s 
displacement in the oxygen polyhedron 28,29 . The proposed system is 
the 0.5PbTiO 3 -0.5(Bi!_ x La x )FeO 3 (abbreviated as the 0.5PT-0.5BL X F, 
0.0 < x < 0.2) whose Bi and La have essentially the same cation 
radius (1.34 A for Bi 3+ and 1.36 A for La 3+ ) however with much 
different ferroelectricity activity. Bi 3+ has the strongest ferroelectric 
activity (0.8A), while La 3+ has the weakest one (0.05A) 28,29 . The aver- 
age bulk CTE can be controlled over a giant range from —7.0 to 
— 0.71 X 10~ 5 °C~\ which contains most known NTE materials. 
The detailed structure was investigated by both neutron powder 
diffraction based on Rietveld refinement and density function theory 
(DFT) calculation. A new concept of spontaneous volume ferroelec- 
trostriction (SVFS) is proposed to elucidate the mechanism of NTE 
control in ferroelectrics. 

Results 

It is known that (l-x)PbTi0 3 -xBiFe0 3 binary system can form a full 
range solid solution with a morphotropic phase boundary (MPB) 
around at x = 0.65, and exhibits unusual physical properties of much 
enhanced c/a, NTE, and Curie temperature (T c ) 23 ' 30 ' 31 . As the x = 0.0 
for the 0.5PT-0.5BL X F, the c/a is significantly enhanced to 1.14 when 
compared with PT (c/a = 1.064). It is the highest in perovskites 
except those metastable compounds needed to be synthesized under 
high temperature and high pressure conditions. The extremely large 
lattice distortion is maintained by the strong Pb/Bi-02 covalency 
which was evidenced by the experimental investigation on the 
electron -charge density distribution 30 . However, such giant lattice 
distortion is significantly reduced by the La substitution (see Supple- 
mentary Fig. S4). With increasing La content, the c axis is much 
reduced, while the a(b) axis is slightly increased. Thus, the c/a is 
significantly reduced. For example, the originally large c/a of 
O.SPT-O.SBL^F is much reduced from 1.14 for the x = 0.0 to 1.057 
for the x = 0.2. Since the coupling of P s and lattice is aligned along the 
c-axis direction, the large decrease in the c axis would imply wea- 
kened ferroelectric property. It is worth noting that even though the 
ionic radius of La 3+ (1.36 A) is essentially the same to Bi 3+ (1.34 A), 
the unit cell volume ( V) apparently contracts by the substitution of 
La 3+ (see Supplementary Fig. S4). The volume contraction is as large 
as 3.1% by the chemical modification from the x = 0.0 to 0.2. Such 
large volume contraction should dominate the NTE property. 

The temperature evolution of unit cell volume for three composi- 
tions of 0.5PT-0.5BL X F (x = 0.0, 0.1, and 0.2) is shown in Fig. 1. It 
is interesting to observe that the NTE is significantly different from 
one another. Without the substitution of La, the 0.5PT-0.5BL X F with 
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Figure 1 | Negative thermal expansion. Unit cell volume of 0.5PbTiO 3 - 
0.5(Bi 1 . x LajFeO 3 (x = 0.0, 0.1, and 0.2) as function of temperature. FE 
and PE mean ferroelectric and paraelectric, respectively. 

x = 0.0 represents a nonlinear and strong NTE in a wide temperature 
range of RT to 600 °C. It is similar to the composition of 0.4PT-0.6BF 
near the MPB 23 . In a lower temperature range (RT ~ 300°C), the unit 
cell volume shows little dependence on temperature with an average 
bulk CTE ofay = -0.83 X 10" 5o C _1 ; However, it contracts dramat- 
ically at elevated temperatures with ot v = -7.0 X 10" 5O C _1 (300 ~ 
600 °C). The average CTE of the overall temperature range (RT to 
600 °C) is ot v = -4.06 X 10" 5 °C _1 . The contraction of unit cell 
volume is as giant as -2.3% for the 0.5PT-0.5BL X F (x = 0.0). It worth 
noting that the NTE of 0.5PT-0.5BL X F (x = 0.0) is among the giant 
NTE materials, such as Mn 3 AN (-1.3%) 17 , BiNi0 3 (-2.5 ~ 
— 3.4%) 9 , and CuO nanoparticles (— 1.1%) U . It also needs to note 
that the present giant NTE in the 0.5PT-0.5BL X F (x = 0.0) occurs in a 
wide temperature range (300 °C) which is broad than other giant 
NTE materials. For example, BiNi0 3 has a colossal NTE over 80 °C 9 ; 
(Mn 0 .96Feo.o4)3(Zn 0 .5Geo.5)N shows a giant NTE over 70 °C 17 ; and 
CuO nanoparticles exhibit a giant NTE over 170 °C n . It is interesting 
to observe that the giant NTE of 0.5PT-0.5BL X F can be distinctly 
weakened by the chemical substitution of La. For the x = 0.1, the 
NTE becomes linear below T c with a CTE of a v = —3.55 X 10" 5O C _1 
(RT ~ 500 °C). With further substitution of La content (x = 0.2), the 
unit cell volume only slightly contracts with increasing temperature. 
A very low thermal expansion is achieved in a wide temperature 
range (RT ~ 400 °C). Here, the CTE is a v = -0.71 X lO^C" 1 
for the x = 0.2. We can achieve that the NTE of 0.5PT-0.5BL X F 
can be tunable over a giant scale of CTE ( — 7.0 ~ —0.71 X 
10~ 5o C _1 ) simply by adjusting La content. It is also important to note 
that such giant tunable range of CTE contains most known NTE 
materials 7 . 

What accounts for such giant-scale control of NTE in the single- 
phase 0.5PT-0.5BL.F? The NTE of the present 0.5PT-0.5BL X F, which 

AV 

is quantitatively described by the CTE in terms of ot v = > is 

AT x Vq 

actually dominated by the volume difference, AV=Vt c — Vrt> 
where Vrt and Vj c are unit cell volumes at RT and T c , respectively. 
As shown in Fig. 1, larger absolute value of AV means stronger NTE, 
while smaller one means weakened NTE. It is known that the unit cell 
volume reaches the minimum value at T c , due to the fact of NTE in 
the ferroelectric phase below T c whereas positive thermal expansion 
(PTE) in the paraelectric one above T c . Here, the Vt c is actually 
dominated by the ionic radius due to the loss of ferroelectricity. 
Above T c , the unit cell volume changes linearly, which is originated 
from lattice thermal vibration. The Vt c of the 0.5PT-0.5BL X F follows 
to a linear extrapolation of unit cell volume in paraelectric phase. 
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However, the Vrt is strong different from one another. One can see 
that the strong NTE of the x = 0.0 is mainly due to considerably 
enhanced V RT , while a near ZTE of the x = 0.2 is due to significantly 
reduced Vrt. The mechanism of NTE control over such giant scale 
would be elucidated according to the detailed studies on chemical 
bonding and hybridization, which will be discussed as follows. 

The crystal structure at RT was refined based on Rietveld method 
by means of neutron powder diffraction. The experimental patterns 
can be well refined according to tetragonal structure model similar to 
PT, and atomic positions could be stably calculated. It is known that 
the equilibrium lattice parameters of ferroelectricity are a balance of 
short-range repulsion favoring cubic paraelectric phase and long- 
range Coulomb force favoring ferroelectric one 32 . P s displacement 
is a key parameter to represent how ferroelectric property changes in 
a unit cell. Figure 2 represents the P s displacements at A-site (Sz A ) 
and 5-site (dz B ) for the 0.5PT-0.5BL X F as function of x. They are both 
reduced at a large content by the substitution of La, indicating a 
weakened ferroelectric property. For example, the dz A is reduced 
from 0.6903 A of the x = 0.0 to 0.4619 A of the x = 0.2 
(Supplementary Table SI). The decrease in the P s displacements 
should originate from the significantly different ferroelectric activity 
of Bi and La atoms. It is known that Bi has the strongest ferroelectric 
activity (0.8 A), whereas La has a very weak one (0.05 A) 29 . The La 
substitution for Bi favors ionic bonding and would weaken the cova- 
lent hybridization of cations with oxygens, which can be also 
reflected from the change in the bond length (see Supplementary 
Fig. S5). With increasing content of La, chemical bonds of A-02 
and B-Ol (the shorter one) are elongated, indicating a weakened 
covalency and hybridization effect. For example, A-02 and B-Ol 
of the x = 0.0 are 2.4287 A and 1.7089 A, respectively. However, 
those of the x = 0.2 are elongated to 2.5103 A and 1.8268 A, respect- 
ively. La substitution favors the formation of cubic paraelectric lat- 
tice, and weakens the coupling between ferroelectricity and crystal 
lattice. Actually, there is a strong correlation between the P s displace- 
ment and the c axis (see Supplementary Fig. S6). The c axis is largely 
reduced by decreasing P s displacement. As a result, the Vr T of 0.5PT- 
0.5BL X F is significantly reduced by the La substitution. Thus, the 
originally giant NTE of 0.5PT-0.5BL X F is strongly weakened by the 
La substitution. 

The effect of La substitution on the covalent hybridization can be 
further supported by density function theory (DFT) calculation. Two 
compositions of 0.5PT-0.5BL X F (x = 0.0 and 0.25) were calculated in 
two 2X2X2 supercells of (Pb 4 Bi4)(Ti 4 Fe4)024 and 
(Pb 4 Bi 3 Lai)(Ti4Fe4)024, respectively. The degree of hybridization can 
be quantitatively described by the value of minimum charge density 
(MED) of specific chemical bonds according to the valence-electron 




Figure 3 | Valence electron-density distributions on the be plane at x = 1/ 
2 of supercell. (a) (Pb 4 Bi4)(Ti4Fe 4 )0 2 4, and (b) (Pb 4 Bi 3 La 1 )(Ti 4 Fe 4 )0 2 4, 
which, respectively, corresponds to x = 0.0, and x = 0.25 of 0.5PbTiO 3 - 
0.5(Bi 1 . x LajFeO 3 through first-principle calculations. The dash line frame 
indicates the original lattice without La substitution. The 0% and 100% in 
color scale corresponds to 0.2 and 8A 3 , respectively, and contours are 
from 0.2 to 4 A" 3 by step of 0.2 A" 3 . 

density distribution (Fig. 3). Without La substitution, there is a strong 
orbital hybridization for the strong NTE composition of 0.5PT- 
0.5BL X F (x = 0.0). The MED value of Bi-02 bond is as high as 0.45 
A" 3 . However, the isovalent substitution of La for Bi weakens the 
hybridization behaviour for all chemical bonds, especially for the 
chemical bond of Bi-02. Here, the MED value is reduced to 0.39 
A" 3 for the composition of x = 0.25 which should have a similar 
low thermal expansion to the x = 0.2. The originally strong ferroelec- 
tric property cannot well be maintained anymore by the incorporation 
of weakly ferroelectric active La for strongly ferroelectric active Bi. And 
thus the P s displacements are all reduced (Supplementary Table S3). 
Therefore, due to the fact of much weakened ferroelectricity and a large 
reduction on the c axis, the Vrt of the 0.5PT-0.5BL X F significantly 
shrinks with increasing content of La. 

Both detailed studies of crystal structure and first principle cal- 
culation have approved that there is a strong correlation between the 
ferroelectric property and the NTE for the 0.5PT-0.5BL X F. The pre- 
sent results further support the conclusion of the correlation between 
ferroelectricity and NTE, which was ever revealed by studies of lattice 
dynamics and temperature dependence of p s 21 > 22 . In order for quant- 
itative characterization on the correlation between ferroelectricity 
and NTE, we introduce a new physical concept of spontaneous 
volume ferroelectrostriction to represent the contribution from 
ferroelectric property to NTE. Spontaneous volume ferroelectrostric- 
tion is defined as co s = — — — x 100%, where V exD and Vn m are, 



respectively, unit cell volumes of experimental and nominal one. The 
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Figure 2 | Spontaneous polarization property, (a) P s displacements at A-site (3z A ) and B-site (3z B ) of 0.5PbTiO 3 -0.5(Bi 1 . x La x )FeO 3 as function of La 
content (x). (b) The sketch of P s displacements of Sz A and Sz B in A-0 12 and B-0 6 polyhedra, respectively. 
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nominal one of V nm is calculated by extrapolation from paraelectric 
to ferroelectric phase. The V nm can be basically considered as the 
volume component contributed from lattice thermal vibration. 
Therefore, the volume difference of V exp — V nni represents the 
quantitative contribution from ferroelectricity. Large value of co s 
means strong NTE, while small one means weak NTE. Similar to 
Landau theory 33 , here we plot the relationship of co s ~ adz 2 A at RT for 
all investigated compositions of 0.5PT-0.5BL X F, where a is the coup- 
ling coefficient between co s and P s (Fig. 4). It is clearly to see that the 
co s shows a good correlation with the square of P s displacement 
(Sz\). On the one hand, strong ferroelectric property induces a large 
value of co s and thus produces a strong NTE such as for the x = 0.0 
without La substitution. On the other hand, weakened ferroelectric 
property induces a small value of co s and produces a low thermal 
expansion such as for the composition of x = 0.2. 

Discussion 

The present study gives a distinct guide to control thermal expansion 
of ferroelectrics. An effective way to control the NTE would benefit 
the CTE matching for different components of devices and compo- 
sites. We can modulate the ferroelectric activity of cations to adjust 
the value of w s and thus control the thermal expansion of ferro- 
electrics. Strong NTE can be achieved by the substitutions of those 
cations with strong ferroelectric activity, such as Bi, Ti, Fe, Zn, and 
Cd, while low thermal expansion by the substitution of cations with 
weak ferroelectric activity, such as La, Mg, Ni, and Zr 28 . The present 
study also offers a useful way to achieve ZTE in the family of ferro- 
electrics. Actually, some of them have been previously experiment- 
ally observed to show such ZTE property, such as (Pb,La)Ti0 3 34 , and 
well-known ferroelectric relaxor of Pb(Mg 1/3 Nb2/3)03 15 . It has been 
known for a long time that perovskite-type Pb-based relaxors exhibit 
near ZTE below the so-called Burns temperature (T B ) 15,35 . Such near 
ZTE property could be understood according to the present concept 
of spontaneous volume ferroelectrostriction. As temperature cooling 
down from T B , the effect of spontaneous volume ferroelectrostriction 
is produced due to the frozen polar nanoregions (PNRs). Therefore, 
the normal thermal expansion is counteracted and thus near ZTE 
appears. Future studies could be focused on the precise control and 
fabrication of ZTE ferroelectric materials for technical application. 
The NTE control could be also extended to other important func- 
tional materials such as multiferroics, since NTE have both strong 
coupling with not only ferroelectric but also magnetic properties. 
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Figure 4 | The correlation between the spontaneous volume 
ferroelectrostriction (co s ) and spontaneous polarization. There is a strong 
correlation of w s ~ ocdz\ for the 0.5PbTiO 3 -0.5(Bi 1 . x La x )FeO 3 
ferroelectrics. 



In summary, NTE of 0.5PT-0.5BLxF single-phase ferroelectrics is 
effectively controlled over a giant range from a strong NTE to a near 
ZTE by the chemical modification. The controllable NTE contains 
most NTE materials. The details of chemical bonds have been studied 
by neutron powder diffraction and DFT calculations. The facile con- 
trol of NTE is strongly correlated to the tunable spontaneous volume 
ferroelectrostriction by modulating the ferroelectric activity of 
cations. Present study establishes an effective way to control thermal 
expansion of functional materials of ferroelectrics. 

Methods 

Sample preparation. Single-phase materials of 0.5PbTiO3-0.5(Bi!. x La x )FeO 3 (x = 
0.0, 0.1, and 0.2) were prepared by the conventional solid-state reaction method. The 
following raw materials were used: PbO, Bi 2 0 3 , Ti0 2 , La 2 0 3 , and Fe 2 0 3 . The well- 
mixed raw materials were calcined at 750° C for 10 h. Then compacted pellets were 
sintered at 1100°C for 2 h covered with the powders of the same composition to 
compensate the loss of PbO and Bi 2 0 3 during high temperature sintering process. 
The surface layers of sintered pellets were removed, and then the sintered pellets were 
carefully crushed into fine powders. 

Structure analysis. Room temperature neutron powder diffraction (NPD) was 
collected at the NIST Center for Neutron Research on the BT-1 high-resolution 
neutron powder diffractometer. The wavelength of synchrotron light was 1.5403 A. 
The initial structural model of PbTi0 3 (space group, P4mm) was adopted for the 
structural refinement performed on software Fullprof. The positions of Pb/Bi/La, Ti/ 
Fe, Ol, and 02 are la site at origin (0,0,0), lb site at (1/2, 1/2, Z B ), lb site at (1/2, 1/2, 
Z 0 i), and 2c sites at (1/2, 0, Z Q2 ) 36 , respectively. The high temperature X-ray 
diffraction (XRD) data were collected on an X-ray powder diffractometer (TTRIII, 
Rigaku, Japan) using Cu Kol radiation. The sample was maintained at a specified 
temperature for 10 min in order to reach thermal equilibrium. The scanning rate of 
angle 26 was 4° min" 1 and the heating rate was 10°C min -1 . The lattice parameters 
were calculated using the software PowderX and TREOR. 

DFT calculation. The electronic structure was calculated based on density functional 
theory (DFT) with generalized approximations (GGAs) performed with Vienna ab 
initio simulation package (VASP), in order to obtain the valence electron- density 
distribution of approximated 2X2X2 supercell of (Pb 4 Bi4)(Ti 4 Fe4)0 2 4 and 
(Pb 4 Bi 3 La 1 )(Ti 4 Fe4)0 2 4, corresponding to tetragonal compositions of 0.5PT-0.5BL X F 
(x = 0.0 and 0.25), respectively. 
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